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ABSTRACT 

The thermal dominant state in black hole binaries (BHBs) is well understood but rarely seen in 
ultraluminous X-ray sources (ULXs). Using simultaneous observations of M82 with Chandra and 
XMM-Newton, we report the first likely identification of the thermal dominant state in a ULX based 
on the disappearance of X-ray oscillations, low timing noise, and a spectrum dominated by multicolor 
disk emission with luminosity varying to the 4th power of the disk temperature. This indicates that 
ULXs are similar to Galactic BHBs. The brightest X-ray spectrum can be fitted with a relativistic 
disk model with either a highly super-Eddington (idisk/^Edd = 160) non-rotating black hole or a close 
to Eddington (Ldisk/^Edd ^ 2) rapidly rotating black hole. The latter interpretation is preferred, due 
to the absence of such highly super-Eddington states in Galactic black holes and active galactic nuclei, 
and suggests that the ULX in M82 contains a black hole of 200-800 solar masses with nearly maximal 
spin. On long timescales, the source normally stays at a relatively low flux level with a regular 62-day 
orbital modulation and occasionally exhibits irregular flaring activity. The thermal dominant states 
are all found during outbursts. 

Subject headings: black hole physics — accretion, accretion disks — X-rays: binaries — X-rays: 
individual (CXOM82 J095550.2-H694047=X41.4-H60) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are non-nuclear 
X-ray sources with apparent luminosities above the Ed- 
dington limit of stellar-mass black holes (BHs). Vari- 
able ULXs are black hole bi naries (BHBs) and may har- 
bor intermediate-mass BHs (IColbert & MushotzkvlilQQl 
iMakishima et al.ll2000[ ). The emission of BHBs has been 
classified into four states based on spectral and timing 
properties: the quiescent, hard, thermal dominant (TP), 
and steep power-law states (McC lintock fc Remillardl 
120061 : [RemiUard fc McClintock 2006i) . The TP state is 
the best understood an d is well described by the st andard 
accretion disk model (jShakura fc SunvaevI I1973D . The 
emergent X-ray spectrum is described by a multicolor 
disk (MCD) model with two parameters: the disk inner 
radius (i?in) and the temperature (Tin) at that radius. In 
the TD state, Rin is constant and the accretion disk is 
thought to extend all the way to the "innermost stable 
circu lar orbit" (ISCO) around the BH (jTakana fc Lewiiil 
|1997[ ). The ISCO radius depends solely on the mass and 
spin of the BH. Therefore, a TD spectrum from an ac- 
creting BH can be used to shed light on its mass and 
spin. 

The TD state has been found during outbursts of many 
BHBs and exhibits specific properties. The disk bolo- 
metric luminosity is I/disk = 47ro'i??jT^* . For constant 
i?in, a 4th power relation between the luminosity and 
the disk inner temperature, Ldisk oc T^, is observed 
(|Kubota fc Makishimall20(M ). Due to up scattering of 
disk photons in the corona, sometimes the observed ex- 
ponent is lower than 4, but can be recovered by applying 
hardening correction (jMcClintock et al.l [20091) . More- 



over, sources in the TD state have low levels of short term 
variability, with very weak or absent narrow band tim- 



broad band power continuum ( 


McClintock & Remillardl 


120061: IRemillard & McClintock! 


20061). 



Spectral surveys of ULXs in nearby g alaxies show 
they are rarely in th e TD state (Feng & Kaaret "200^ 
Stobbart etall 120061 : [Winter .etjjL 2006 : Sor ia et afl 



20091 ) ■ The spectra of a few ULXs are super-soft and can 



be modeled by an MCD with an inner temperature close 
to 0.1 keV. However, repeated observations revealed an 
evolution pattern inconsistent with the Ldisk oc rela- 
tion rulin g out inter pretation as thermal dominant disk 
emission (iLiul I2008D . Some other ULXs show soft ex- 
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cesses in the spectra that could be fitted by c ool disk 
emission (jKaaret et all 120031: iMiller et all I2004D . How- 
ever, the cool disk is not the dominant component in 
their spectra, thus, the sources are not in the TD state. 
Finding a ULX in the TD state would demonstrate that 
the ULX has emission states similar to Galactic BHBs 
and allow inference of the BH mass and spin. 

The starburst galaxy M82 contains one of the most lu- 
minous ULXs CXOM82 J095550-h694047 ( =X41.4-f60; 
iKaaret et all 120011: iMatsumoto et al.ll200 1'). in nearby 
galaxies. The ULX was first identified with Chandra 
at a lumino sity higher than the Eddington limit of a 
5OOM0 BH (|Kaaret et al.|[200ll:[MaFsumoto et al.|[200l . 
On the sky, it lies near a super star cluster, which 
was speculate d to be the b irth place of a massive BH 
{Portegies Zwa rt et al.ll2004 ). Low frequency QPOs and 
broadband timi ng noise, detected in th e central region of 
M82 (Strohmaver fc Mushotzkvl 120031 : iDewangan et all 
l200llMucciarelli et al.ll2006ir and later confirmed to orig- 
inate from this ULX ( Feng fc Kaaretll 2007f). sug gest that 
the ULX harbors a massive BH. King_&Delmen (2005|) 
also suggests that the ULX contains a massive BH which 
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is the nucleus of a satellite galaxy merging with M82. 

Positive identification of the emission states requires 
both timing and spectral information. Here, we describe 
simultaneous observations exploiting the high angular 
resolution of Chandra to isolate the ULX spectrum from 
diffuse emission and nearby sources and the large collect- 
ing area of XMM-Newton to obtain timing information. 
We use the joint timing and spectral information to iden- 
tify the ULX emission states and unveil the nature of the 
source. 

2. RESULTS FROM JOINT CHANDRA/XMM-NEWTON 
OBSERVATIONS 

Three Chandra (ObsIDs 10027, 10025, and 10026) and 
XMM-Newton (ObsIDs 0560590101, 0560590201, and 
0560590301) observations were performed on 2008 Oc- 
tober 4, 2009 April 17 and 29, respectively. Each Chan- 
dra observation has an effective exposure of about 18 ks 
and overlapped with an XMM-Newton observation that 
lasted 30, 40, and 50 ks, respectively. 

To reduce the effects of pileup, the Chandra observa- 
tions placed the ULX 3.57" off the optical axis, so the 
source covers multiple pixels, and the CCD was operated 
in a l/8th sub-array mode to reduce the readout frame 
time. Using the pileupjnap tool in CIAO 4.1.2, we mea- 
sured that the ULX resulted in 0.25 counts per frame 
(~13% pileup fraction) on the brightest 3-by-3 pixel is- 
land in the observation on 2008 Oct 4, and found that 
about 85% of the total events were recorded on pixel is- 
lands with counts per frame larger than 0.06 (^3% pileup 
fraction). Therefore, pileup must be taken into account 
in the spectral fitting. We note that the other two obser- 
vations have similar, mild pileup, which can be corrected 
by applying a pileup model in the spectral fitting. 

Each Chandra spectrum was extracted from a source 
region found by wavelet detection. Instrument responses 
were calculated using the calibration files in CALDB 
4.1.3. The spectrum channels were re-grouped by a fac- 
tor of 8 in 0.3-1 keV, 4 in 1-4 keV, 8 in 4-6 keV and 16 
in 6-8 keV, respectively and were fitted in XSPEC 12.5 
to a power-law model or MCD model (DISKBB) sub- 
ject to interstellar absorption and CCD pileup with back- 
ground subtracted from a nearby, source free region. To 
minimize contamination from diffuse emission and piled 
events, we fitted the data in the energy range from 0.7 
to 7 keV. During the fits, the pileup grade migration pa- 
rameter a was free, and /psf (the fractional events in the 
spectrum being treated for pileup) was frozen at 0.85 
based on the estimate discussed above and is also the 
best-fit value found from the fits. 

We found that the MCD model provides an adequate 
fit to the data, while the power-law model does not. 
No reasonable pileup parameters can be found with the 
power-law model; the a parameter always goes to zero, 
indicating the pileup fraction in the spectrum is zero, 
which is contradicted by the pile-up estimates. For the 
2008 Oct 04 observation, the fit is not acceptable with 
= 102.1 for 70 degrees of freedom (dof). If we fix a at 
its best-fit value found from the disk model, the goodness 
of fit is even worse, with x^/dof = 201.8/71. In contrast, 
the MCD model provides a good fit xVdof = 67.8/70 
and a reasonable value of a. Figure [T] shows the spec- 
trum from the Chandra observation on 2008 Oct 04 and 
the comparison between the MCD and power-law model. 



Table 1 

Spectral parameters of X41. 4-1-60 fitted with the MCD model 
subject to interstellar absorption and CCD pileup. 



Parameters 


2008 Oct 04 


2009 Apr 17 


2009 Apr 29 


a 

Nil 

Tin 
Rin 

fx 

-t'disk 


0.18 ± 0.05 
1.13 ±0.04 
1.52 ± 0.09 
339 ± 24 
5.3 ±0.8 
7.9 ±2.0 
67.8/70 


0.26 ±0.06 
1.07 ±0.04 
1.44 ± 0.09 
330 ± 28 
4.3 ±0.5 
6.0 ± 1.8 
79.0/68 


> 0.29 
1.05 ±0.07 
1 10+"-^^ 

011+24 
1 9+0.27 

^•'^-o.os 

1 0+1.4 
1-8-1.0 

42.5/53 



Note. — a is the pileup grade migration probability. A^h is 
the absorption column density in 10^^ cm~^. 7]^ is the disk inner 
temperature in keV. Hin is the disk inner radius in km calculated 
assuming a distance of 3.63 Mpc and a disk inclination angle of 
zero, /x is the absorbed flux in 1-7 keV corrected for pileup in 
10~^^ ergs~^ cm~^ s~^. I/disk is the bolometric disk luminosity 
in 10*" ergs s~^. All errors are quoted at 90% confidence level. 
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Figure 1. The Chandra data spectrum of the ULX on 2008 Oct 04 
and comparisons between the MCD and power-law model, (a) The 
data spectrum and the best-fit MCD model modified by interstellar 
absorption and CCD pileup (see Table[T]for parameters), (b) Data 
to model ratio for the MCD model (x^/dof = 67.8/70) (c) Data 
to model ratio for the power-law model {x^ /dof = 201.8/71) with 
the pileup parameter a fixed at its best value found from the MCD 
model, (d) Data to model ratio for the power-law model (x-^/dof = 
102.1/70) with a free during the fit, which, however, converged to 
an unphysical value [a = 0) indicative of a zero pileup fraction. 
Therefore, the MCD model is favored and the power-law model is 
ruled out from the spectral fitting. 
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Figure 2. Bolomctric disk luminosity versus the disk tempera- 
ture at inner radius of X41.4+60 derived from the three Chandra 
observations by fitting the data to the MCD model subjeet to in- 
terstellar absorption and CCD pileup. A constant disk inner radius 
is derived from the fits, as a consequence of which, the disk lumi- 
nosity is scaled with the inner temperature to a power of 4. The 
dashed line indicates a best-fit 4th power-law relation to the data. 
The errors are at the 1 a level. 



Please note that the data/model ratios from the two 
power-law models are not evenly distributed. 

The MCD model also provides better fits than power- 
law for the other two observations, see Tabled] For the 
observation on 2009 Apr 17, 26.3 of arise from in- 
strumental response features on individual channels at 
2.1 and 3.2 keV. Ignoring these two channels does not 
change the spectral parameters at all. We excluded them 
because we are only interested in the broad continuum. 
For the observation on 2009 Apr 29, possible emission 
line features are seen at energies from 4.5 to 7 keV. For 
the same reason, we only use 0.7-4.5 keV for this obser- 
vation. We note that for the other two observations the 
spectral parameters derived from the 0.7-4.5 keV range 
are completely consistent with from the 0.7-7 keV range. 
Although adding a power-law component to the model 
does not improve the fits, we did that to place an upper 
limit on the power-law flux with the power-law photon 
index bounded in the range 2.1-4.8 as found in the TD 
state. The upper limits on a power-law component in 
2-20 keV are estimated to be 4%, 17%, and 20%, respec- 
tively. This is consistent with the definition of the TD 
state. 

Power spectra were calculated from XMM-Newton 
data using combined PN and MOS events in a circular 
region around the source with a radius of 18" from the 
common good time intervals without timing gaps (foUow- 
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Figure 3. Lightcurve of the whole M82 galaxy observed with 
RXTE from 2008 Feb 5 to 2009 Aug 2. The crosses are the observed 
X-ray flux in the energy range of 2-10 keV, which is dominated 
by X41.4+60. The source spends most of its time at a relatively 
low flux level, around which the flux is modulated at a period of 
62 days, and occasionally shows flaring activity lasting days to 
months. The three arrows indicate the times of the joint Chandra 
and XMM-Newton observations, in which QPOs disappeared and 
the source was identified as in the TD state. 



ing lFeng fc KaareO 120071) . We searched for QPOs from 
1 to 1000 mHz in various energy ranges, but did not 
find significant signals. We repeated the search in the 
half of the circular source region where contamination 
from source s other than the ULX is minimized (region A 
in Fig. 2 of lFeng fc Kaaretll2007[) . but still found noth- 
ing. The power spectra are consistent with that from the 
white noise. The 3(7 upper limits of the total noise power 
(in units of rms/mean) in 1-1000 mHz in the energy range 
of 0.3-10 keV are 6.3%, 6.6%, and 6.9%, respectively from 
the three observations. 

These XMM-Newton observations are able to detect 
QPOs at the strengths previously observed. The disap- 
pearance of QPOs and the low noise level suggest that 
the source was not in the hard state. The Chandra en- 
ergy spectra show that they are incompatible with a 
power-law form, but can be adequately modeled with 
a dominated MCD model with a constant inner radius. 
Figure [2] shows that the bolomctric luminosity from the 
accretion disk varies with the 4th power of the tem- 
perature. The correlation coefficient between log idisk 
and log Tin is 0.9995 with a chance probability of 0.02 
from the three points. All of the emission properties, 
low timing noise, spectrum dominated by an MCD with 
low power-law fraction, and luminosity proportional to 
disk temperature to the 4th power are consistent with 
the source being in the TD state and are strong evidence 
in favor of the interpretation of the energy spectrum as 
emission from an optically thick accretion disk. 

The lightcurve from RXTE for the whole M82 galaxy 
is shown in Figure [3] from 2008 February 5 to 2009 Au- 
gust 2. The Chandra/XMM-Newton described here all 
occurred when the fiux from M82 was well above its nom- 
inal level. The ULX was, by far, the brightest source in 
the three Chandra observations. Therefore, these out- 
bursts from M82 are dominated by the ULX, with a mi- 
nor contribution from another source (X42. 3-1-59) active 
only in the first joint observations. Thus, the TD states 
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were all observed during bright outbursts of the ULX. 

Only two other Chandra observations have manage- 
able pileup. An observation on 2005 February 4 revealed 
a relatively low X-ray flux and a hard pow er-law spec- 
trum inconsistent with an MCD model p<!aaret et alj 
|2006[) suggesting the hard state. We previously analyzed 
an observation from 2007 July 2, but fixed /psf = 0.38 
based on the assumption that pileup was important only 
in the 3-by-3 pixel cluster with the highest count rate 
(jKaaret et al.|[2009l ). Utilizing the new pileup jnap tool, 
we realized that /psf was underestimated. Using the im- 
proved procedures described above, we find xl — 2.2 
for a power-law model and xt = 1-4 for a MCD model 
with parameters similar to those from 2008 Oct 04. Lack 
of timing information precludes firm conclusion, but the 
TD state is suggested. 

2.1. Fitting with a Relativistic Disk Model 

The standard MCD model does not include relativis- 
tic effects. Such effects are particularly important for 
rapidly spinning BHs. For accurate estimation of the BH 
properties, we fitted the spectrum from when the source 
was brightest, the 2008 obser vation, with a fully relativis- 
tic MC D model KERRBB (iLi et al.ll20"05l). W e set zero 
torque ()Li et al.ll2005t iMcClintock et al.ll2006D at the in- 
ner edge of th e disk, the distance to the source as 3.63 
Mpc jF rcedm an et al.|[T9 94'). the hardening factor as 1.7 
(jMakishima et al.ll2000D ~ and switched on self-irradiation 
and limb- darkening. The pileup parameters were fixed at 
their best values found from the MCD model. This left 
five free parameters in the model: the disk inclination, 
BH mass, spin, accretion rate, and the interstellar ab- 
sorption column density. 

To test if the BH is maximally spinning, we fixed the 
spin parameter at two extreme values a* = and 0.9986, 
respective ly, representin g non-spinning and maximally 
spinning (|Li et all [2001 . For a non-rotating BH, the 
accretion rate exceeds the Eddington limit by a factor 
of 160. For Galactic BHBs, the observed Eddington ra- 
tio is as high as 7 in V4641 Sgr despite the unc ertainty 
on the distance (jMcClintock fc RemillardI 120061 ). When 
the luminosity is high and the radiation pressure is domi- 
nant, the thin accretion disk is predicted to be able to ex- 
ceed the Eddington limit by a factor of a few (jBegelmanl 
12002^). Therefore, we adopt an upper limit of 10 for the 
Eddington ratio as being physically reasonable. We note 
that active galactic nuclei do not appear to exceed the 
Eddington limit by more than a factor of 10 in any accre- 
tion state (jWarner et al.|[200l) . The Eddington ratio for 
a non-rotating BH is too large and this scenario should be 
ruled out. In contrast, if the BH is maximally spinning, 
the derived Eddington ratio is only 2. 

Freezing the spin at the maximum allowed value re- 
sulted in a best-fit BH mass of 660 solar masses and a 
90% error range of 300-1250 solar masses, a disk inclina- 
tion angle of at least 60°, and xVdof = 67.7/70. The 
other two observations produce similar results. We also 
applied the model to the three observations simultane- 
ously, imposing the same mass, spin, and inclination, and 
allowing the mass accretion rate and absorption column 
density to vary individually. Fixing the spin parameter 
a* = 0.9986 led to a best-fit BH mass of 500 solar masses, 
a 90% error range of 300-810 solar masses, an inclination 
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Figure 4. Energy spectra of X41. 4-1-60 from the three Chandra 
observations and the simultaneous fit to the KERRBB model sub- 
ject to interstellar absorption and CCD pileup. The spin parameter 
is fixed at its maximum a* = 0.9986. For other parameters please 
refer to the discussion of this model in the text. The residues are 
for the observations on 2008 Oct 4, 2009 Apr 17, and 2009 Apr 29, 
respectively, from top to bottom and in units of a. 

of 59-79 degrees, and x^/dof = 189.0/195; see Figure H 
The Eddington accretion rate at the highest luminosity 
observed is 2.5. The inferred Eddington ratio exceeds 10 
for a* < 0.93, which we took as the lower limit of the 
spin. Fixing a* = 0.93, the best-fit BH mass is 430 solar 
masses with a 90% error range of 190-570 solar masses, 
and the inclination is larger than 60°. Thus, the spectral 
fitting suggests that the ULX contains a rapidly spin- 
ning (a* > 0.93) intermediate mass BH of 200-800 solar 
masses. 

3. DISCUSSION 

Joint Chandra and XMM-Newton observations have 
enabled us to precisely measure the spectral and timing 
behavior of this ULX in M82. Disappearance of QPOs 
and low timing noise, an X-ray spectrum best-fitted with 
an MCD model, and an idisk oc T-^^ pattern in the spec- 
tral evolution provide strong motivation to interpret this 
behavior in terms of the TD state seen in stellar-mass 
BHBs. The MCD model indicates that the source spec- 
trum shows spectral curvature at a few keV and thus 
could also be fitted by more complicated models like 
a cool optically thick corona which has been proposed 
recently as an indicato r of a new 'ultraluminous state' 
(iCladstone et al.ll2009| ). However, a transition from the 
hard spectrum (incompatible with the steep power-law 
state) seen at lower luminosities to the ultraluminous 
state would be difficult to explain and we, therefore, do 
not discuss this proposed state in detail. 

Simply scaling the Ldisk oc pattern to those from 
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Galactic BHBs suggests that the source contains a more 
massive BH, because the compact object mass is pro- 
portional to the square root of the disk luminosity if 
two BHs have the same disk inner temperature and spin 
(jMakishima et al1l2000[ ). Fitting with a fully relativistic 
MOD model leads to a consistent result. A conservative 
estimate of the BH mass by allowing an Eddington ra- 
tio of as high as 10 suggests a BH of at least 200 solar 
masses but less than 800 solar masses. This is coincident 
with the theoretical estimate of the mass of this BH if it 
is created in a nearby star cluster by runaway collisions 
(jPortegies Zwart et al.ll2004l) . The fast spin of the source 
makes it efficient in extracting gravitational energy. 

The spectral fitting also indicates a relatively high 
inclination angle, 59-79 degrees, of the accretion disk. 
For a maximally spinning BH of a few hundred solar 
masses, the relativistic effect has largely eliminated the 
limb-darkening effect; viewing the disk at a high incli- 
nation angle receives almost the same flux as at a low 
angle. Therefore, the observed high flux is not a prob- 
lem for a nearly edge-on disk. The X-ray flux from M82 
is modulated at a period of 62 days, interpreted as due 
to orbital motion of this ULX binary, which must con- 
tain a giant or super-giant star inf erred from the pe- 
riod assum ing Roche-lobe overflow (|Kaaret et al .1120061 : 
iKaaret fc' Fcng 2007). With knowledge of the BH mass 
(assuming 200-800 solar masses) and the 62-day binary 
orbital period, the binary separation of the ULX can be 
calculated with the assumption of Roche lobe overflow 
as (3 — 4) X 10^^ cm, which is insensitive to the mass of 
the companion star. This is larger than all separations 
of low-mass BHBs with a dynamical measurement of the 
mass, and larger than that of GRS 1915+105, which has 
the largest separation known so far, by a factor of a 
few. This ULX and GRS 1915+105 share a similarity 
that they have both been active for many years with- 
out quenching to the quiescent state. The large separa- 
tion and consequently a huge reservoir of accretion mass 
could be the factor that determines their long-term ac- 
tivity (Dccgan ct al. 2009). 
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and the Foundation for the Author of National Excellent 
Doctoral Dissertation of China under grant 200935. PK 
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